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ABSTRACT 

We investigate the possibility of using cosmological observations to probe and constrain 
an imperfect dark energy fluid. We consider a general parameterization of the dark 
energy component accounting for an equation of state, speed of sound and viscosity. We 
use present and future data from the cosmic microwave background radiation (CMB), 
large scale structures and supernovae type la. We find that both the speed of sound and 
viscosity parameters arc difficult to nail down with the present cosmological data. Also, 
we argue that it will be hard to improve the constraints significantly with future CMB 
data sets. The implication is that a perfect fluid description might ultimately turn out 
to be a phenomenologically sufficient description of all the observational consequences 
of dark energy. The fundamental lesson is however that even then one cannot exclude, 
by appealing to observational evidence alone, the possibility of impcrfcctncss in dark 
energy. 



1 INTRODUCTION 



, The observed acceler ated expansion of the universe (see 
e.g. l|Astier et alj|2006h ) is usually ascribed to the existence 
of a cosmic fluid with a negative pressure comparable to its 
small but cosmologically significant energy density. While 
numerous po ssible origins have bee n proposed for this dark 
energy (see (I Cop eland et alj l2006t ) for a review) most of 
them are variations of the cosmolog ical constant or the 
scalar field scenario (|Wetterichl [l988). These include sev- 
eral different kinds of c ouplings to matt e r (baryons, neu- 
trino s , cold dark matter jAmendolal [20o5; iBrookfield et al.l 
l200d : iMota fc Shawl [20071 . l200fj ; lKoivistoll2005l )), couplings 
to gravity (the curvatur e scalar, Lovelock invariants, etc 
l|Koivisto fc Motall2007bl lcl)). and also non-canonical kinetic 
terms (phantoms, tachyons, K -essence, etc (|Gibbonsll2002l ; 
lArmendariz-Picon et alj|200ll )). With such a zoo of models, 
it is important to investigate, and search for, some particu- 
lar feature or property of dark energy, which could be used 
to rule out some of (or at least distinguish among) all these 
possibilities. 

It is already well known in the literature that to discrim- 
inate between different classes of models, it is not sufficient 
to consider just the background expansion. One has also to 
study the evolution of cosmological perturbations. General 
Relativity dictates that any component with its equation of 
state w 7^ — 1, where w — — 1 corresponds to the cosmo- 
logical constant, must fluctuate. Therefore a generic dark 
energy component has perturbations which couple to mat- 
ter perturbations. These however can be small, since if one 
has w « —1, the component can be nearly smooth, and if the 
Jeans length of dark energy is large, its perturbations may be 
confined to very large scales only. This is typical for the usual 



minimally coupled quintessence models, since their sound 
speed of perturbations, cf am , is equal to th e light speed, 
cf„ m = 1, which set s a large Jeans length (|Bean fc Dord 
l2004l ; lxia et al.ll2007l ). However , for some other dark energy 
candidates that is not the case (IMota fc van de Bruckll2004l ; 
iBagla et all2003l ; iBento et al.ll2002l ). and such feature might 
help to differentiate between these classes of models. 

In addition to the w and c^ am , there is an impor- 
tant characteristic o f a gener al cosmic fluid which is its 
anisotropic stress a (H3I1998). This vanishes for a mini- 
mally coupled scalar field and perfect fluid candidates, but 
is a generic proper ty of realistic fluids with finite shear 
viscous coefficients l|Schimd et al" l l200rj ; lBrevik et al. 1 12004 
iNoiiri fc Odfnt sov 2005). Basically, while w and cf am deter- 
mine respectively the background and perturbative pressure 
of the fluid that is rotationally invariant, a quantifies how 
much the pressure of the fluid varies with direction. In fact 
the anisotropic stress perturbation is crucial to the under- 
standing of evolution of i nhomogeneities in the early ; radia - 
tion dominated universe (|HvJll998l : Ikoivisto fc Mot'all2006h . 
Therefore an obviously interesting question is whether 
present observational data could allow for an anisotropic 
stress perturbation in the late universe which is dominated 
by the mysterious dark energy fluid l|lchiki fc Takahashil 
l2007tlBajrowlll997h . 

The cosmological effects of a due to possible viscosity 
of dark energy are however quite neglected in the literature. 
The main reason for disregarding the anisotropic stress in 
the dark energy fluid might be that conventional dark en- 
ergy candidates, such as the cosmological constant or canon- 
ical scalar fields, are perfect fluids with a = 0. However, 
since there is no fundamental theoretical model to describe 
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dark energy, there are no strong reasons to stick to such as- 
sumption. Moreover, couple d scalar fields have in deed a non- 
negligible anisotropic stress l)Schimd et al.l [2006). and dark 
energ y vector field candidat e s (which have been proposed 
in (lArmendariz-Piconl |2004| ; iKiselevI 12004 IZimdahl et"aH 
120011 )1 also have a ^ 0. 

The aim of the present paper is to investigate the po- 
tential cosmological signatures of a very general dark energy 
component, which is characterize by an equation of state, 
sound speed and anisotropic stress. In section [2] we review 
the parameterization of a generalized cosmological fluid and 
comment its relation to some recent studies of anisotropies 
dark energy. The parameterization will then be subjected to 
the most detailed and most extensive scrutiny this far. The 
data and method utilized for this are described in section [3] 
In the section [4] we use the most recent cosmological data 
to constrain the properties of dark energy. Section [5] is de- 
voted to investigate how much the future data could be able 
to improve the constraints. We conclude by stating the fun- 
damental uncertainty in the properties of dark energy but 
also mention some cases where a positive detection could be 
established. 



2 DARK ENERGY STRESS 
PARAMETERIZATION 

Consider a general fluid with the energy momentum tensor 



Tpv = pu^Uv + ph^v + E,, 



(1) 



where is the four-velocity of the fluid, and the projection 
tensor /i M „ is defined as h^v = g M „ + u^Uv. Here E M „ can 
include only spatial inhomogeneity. At the background level, 
the evolution of the fluid is determined by the continuity 
equation, 



p + 3H(l+w)p = 0. 



(2) 



The effects to the overall expansion are therefore determined 
by the equation of state w alone. We define a perfect fluid 
by the condition £ M „ = 0. The condition for the adiabaticity 
of a fluid is p = p(p) , which implies that the evolution of the 
sound speed is determined by the equation of state alone. 
Generally, however, the sound speed is denned as the ratio 
of pressure and density pert urbations in the frame comoving 
with the dark energy fluid (|Weller fc Lewisll2003h , 

Sp 



2 



dp 



\de 



(3) 



In the adiabatic situation one has c 2 a 



dp I dp 



w — 3H ^ +ul j , but in general the sound speed is an indepen- 
dent variable. In the following we will consider a constant 
equation of state for simplicity. 

Taking these considerations into account, the evo- 
lution equations for the dark energy density perturba- 
tion S and velocity poten tial 9 in the synchronous gauge 
l|Ma fc Bertscmngenll995h . can be written as 

6 = -(l+w)^[k 2 + 9H 2 (c 2 am -w)] ^ + ^| 



3-W(cf am - w)8 
-H(l-3cf am )e + 



2 7,2 
1+W 



(4) 
(5) 



where h is the trace of the synchronous metric perturbation. 
Here a is the anisotropic stress of dark energy, related to no- 
tation of Eq.JTJ by (p + p)a = — (kikj — |<Sij)£ lJ . From the 
above equations it is then clear that, while w and cf am deter- 
mine respectively the background and perturbative pressure 
of the fluid that is rotationally invariant, a quantifies how 
much the pressure of the fluid varies with direction. 

To close the system of equations, we describe the evolu- 
tion ofthe_anisotropic stress with an equation adopted from 
Hu ([Hulll998T ). 



■ i on 8 c; 
a + SHa = — 



3 1 + w 



h 

2+3.). 



(6) 



This parameterization leads to reasonable results and ap- 
proximates the evolution of any fluid present in the standard 
cosmological model, in particular neutri nos and photons 
which have a non-zero anisotropic stress l|Hu fc Eisensteinl 
1999). More specifically, for those relativistic components, 
the correct choice for the viscous parameter is c vis = 1/3. A 
perfect fluid (vanishing shear viscosity) should have c 2 is = 0. 

Equivalently, one can describe the physical properties of 
a by introducing the rescaled parameter a V i S = c 2 ia /(l + w). 
While c 2 ia is somehow physically analogous to sound speed 
squared, the a V i a is the quantity directly multiplying the 
source term of the stress (see RHS of eq. ©). In this study 
we will use both the c 2 ia and a v i s parameterizations. Since 
w is constrained near w = —1, where the relation between 
these parameters is divergent, using one or the other might 
lead to different results and interpretations. The statistical 
details also depend on which parameter one assumes a uni- 
form distribution, and it is useful to test how robust ones 
conclusions are to such assumptions. 

Usually, a dark energy fluid with nonzero a V i S generates 
shear stress which tends to smoothen its distribution. How- 
ever, the consequences to phantom dark energy are qual- 
itatively different and for such a fluid, with w < —1, a 
shear stress drives the clustering. With negative a v i 3 , ex- 
ponential growth is typical for all kinds of dark energy. In 
the case of positive a V i S , and constant w and cf arn , the ef- 
fects are confined to sup erhorizon scales and typically small 
l|Koivisto fc Motal [2006 ) . This is in contrast to anisotropic 
stress which originates from modifications of the gravity sec- 
tor or quintessence couplings, since they typically modify the 
gravitational pot entials at small scales . A recent study by 
Amendola et al l|Amendola et al.1 [2004) considers the pos- 
sibility of using weak lensing to obtain limits on the dark 
ene rgy parameters motiv ated partially by modified grav- 
ity (I Amendola et al.1l2004r i. Our study can be considered as 
an exploration of complementary aspects since the effects 
we encounter, occur (except for negative cf am or a V i S ) at 
some orders of magnitude larger scales than those possible 
to probe with weak lensing experimen ts. An interesting ap - 
proach is also that of Caldwell et al ([Caldwell et al.|[2003h 
who, parameterizing directly the deviation from the general 
relativistic perfect fluid metric, and assuming it to depend 
on the amount of dark energy, find that this assumption 
leads to relatively weak constraints on cosmological scales, 
whereas the bounds on the solar system are known to be 
tight. 

To illustrate the effect of c 2 is on the CMB power spec- 
trum, we have in Figure [1] shown the CMB temperature 
power spectra for models with different values of w and 
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Figure 1. The CMB temperature power spectrum for models 
with different values of w and c 2 is . The upper panel is for a 
model with w = —1.2 and c 2 . 0, while the lower panel shows 
the results for a model with w = —0.8 and c 2 is 0. The grey 
shading indicates the cosmic variance around the models with 
c vis = 0' A" models plotted here have c 2 am = 1. 

c 2 i3 . In models with w = —1.2, one sees that the deviation 
from the <? viB — model starts to become significant with 
c vis — —25. With w — —0.8, the deviations remain small, 
even with <? via = 1000. From this we would expect that it 
is possible to find some lower limit on c 2 is in models with 
w < — 1, while it will be difficult to find upper limits on (? via 
in models with w > — 1. 

In Figure we have plotted the CMB temperature 
power spectra for models with different values of w and cf am . 
We see that for values of cf am increasing over cf am = 1, the 
shape of the power spectra seems to stabilize, even when 
cf am S e t s as large as cf am = 100. So also in the case of 
cfam > it looks difficult to distinguish between the differ- 
ent models. Such difficulty can be understood noticing that, 
increasing cf am will lead to a smoothening of a quintessence- 
like case. This is itself already almost smooth. Therefore we 
do not see much difference even when cf am = 100. Notice, 
however, that negative cf am is also conceivable, and it corre- 
sponds to gravitational instability. This leads to the cluster- 



Figure 2. The CMB temperature power spectrum for models 
with different values of w and c 2 am . The upper panel shows the 
results for a model with w = —1.2, while the lower panel shows 
models with w = —0.8. The grey shading indicates the cosmic 
variance. All models shown here c 2 . =0. 



ing of dark energy and thus to possibly observable effects. 
In Figure [2] we see that the observable effect of cf am < is 
much more evident than for c 2 am > 0. Such possibility will 
also be investigated in the following sections. 



3 DATA AND METHODS 

In our analysis we use data from observations of the 
anisotropies of the Cosmic Microwave Background, the 
distribution of large scale structures (LSS), the distance- 
redshift relation from type la supernovae (SNIa), and an 
additional prior on the Hubble parameter. The parameter 
estimation analysis is done using a modifie d version of the 
Mark ov chain Monte Carlo code CosmoMC l|Lewis fc Bridle! 
120021) . 

In our basic cosmological model, we assume a flat, ho- 
mogeneous and isotropic background spacetime, and a sim- 
ple power law primordial power spectrum. Thus, we allow 
the six parameters {Qth 2 , Q m , log^O 10 ^), h, n B ,r} to vary 
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freely. Here Qi is the ratio of the energy component i to the 
critical density today, As is the amplitude of the primordial 
power spectrum, h is the dimensionless Hubble parameter 
today, n 3 is the scalar spectral index and r is the optical 
depth at recombination. The exact definitions of the param- 
eters are given by the CosmoMC code. In addition to these 
basic six parameters we have also studied changes in w, e 2 is 
(or a vis ) and cf am . 

The most constraining CMB dataset at present is the 
3-year data release from the WMAP team. In our analysis 
we have used bot h the temperatur e dHinshaw et al.l [20071 ) 
and polarization (|Page et al.l l2007h data from this exper- 
iment together with the likelihood code provided by the 
WMAP team 0. Since the effects we are studying here are 
most prominent on large angular scales in the CMB signal, 
we do not take additional small-scale CMB experiments into 
account in this analysis. 

In one case we have also utilized LSS data from the 
Sloan Digital Sky Survey Lumino us Red Galaxy Sample 
(SDSS-LRG) (jTegmark et alj l2006h . When using informa- 
tion on SNIa distance-redshift relatio n, we use data fro m 
the Supernova Legacy Survey (SNLS) (lAstier et alj|2006h . 

We have also used a prior of h — 0.72 ± 0.08 
from the Hubble Sp ace Telescope Key Project (HST) 
l|Freedman et al.ll200ll ) and a top- hat prior on the age of the 
universe, 10Gyr<Age<20Gyr, throughout the entire analy- 
sis. 



4 CONSTRAINTS FROM PRESENT DATA 
4.1 Constraining dark energy viscosity 



To start with we have used the c v 
ing at two distinct scenarios. In 



, s parameterization, look- 
one case we have {w > 



-l,c vis > 0}, and in the other case {w < — l,c vis < 0}. 



Here, and for the rest of subsection l4.ll we have set c larn = 1. 

In Figure [3] we show the 65% and 98% confidence level 
(CL) contours in the uj— c 2 is plane for the two main scenarios 
mentioned above. We see that in the first case there is no 
clear degeneracy between c 2 i3 and w, and we find no upper 
limit on c 2 is here. This is exactly what we would expect 



-1 c 2 - 

"-"in. 



> o}, 



from an inspection of Figure [T] when {w > 
the effects of varying c 2 is are negligible. 

For the latter case, where w < — 1 and c 2 is < 0, we do 
however find a degeneracy between these parameters, and 
in this case we also find lower limits on c 2 is . Using only 
WMAP data we find c 2 via > -19.5 at 95% confidence level 
(CL) from the 1 dimensional distribution of c 2 is . Since there 
are observable effects of a non-zero c 2 is in this scenario, we 
also tried to add SDSS-LRG and SNLS data to see if this 
would improve our constraints on c 2 is . In this case the lower 
limit on c 2 is becomes c 2 is > —24.9 at 95% CL. We can 
also easily see from the figure that adding SDSS-LRG and 
SNLS data does not improve the limits on c 2 i3 , only on w. 
The reason is that the effects of c 2 is occur at larger scales 
than probed by the SDSS-LRG. Therefore this additional 
data would only be interesting if it could break degeneracies 
between c 2 is and other parameters that govern the shape of 
the power spectrum for low Is. We see from Figure [T] that 



http://lambda.gsfc.nasa.gov; version v2p2. 



-0.4 
-0.5 
-0.6 
-0.7 
-0.8 
-0.9 



-1 







95% 















50 



,100 



150 



200 



-1 



-1.1 



-1.2 



-1.3 



-1.4 



-1.5 



\ \^ 68% 




\. 95% 








\ 95% 







-25 



-20 



-15 



-10 



Figure 3. 68% and 95% contours of c 2 . and w. In the upper 
panel w > —1 and c 2 > 0, while in the lower panel w < —1 and 
c vis ^ W e see that the constraints are much tighter for negative 
c 2 is than for positive. In the right plot we have also shown how 
the constraints improve when including SDSS-LRG and SNLS 
data (dotted black contours). In these models c 2 am = 1. 



although we do have a degeneracy between c 2 is and w, the 
limits on c 2 is would only be improved by some data set that 
would favor a w < — 1. 

Next, we turn to the a V i S parameterization. To con- 
strain this parameter we focus on a model where a V i S < 0, 
since we find no upper limits on a positive a V i S . In Figure 
[4] we show the 65% and 98% CL contours for w and a V i S 
in this model. We see that a V i a is bounded from below and 
that values of a v i s slightly below zero are allowed. In this 
case we find that a V i S > —0.23 at 95% CL. 

Since we got a lower limit on c 2 is when having w < — 1, 
we would also expect to find an upper limit on a V i S in this 
case. In Figure [5] we show the 1 dimensional marginalized 
likelihood for log a V i S for a model with w = —1.2. We find 
that loga„ is < 3.2 at 95% CL (a v i s < 24.5), which corre- 
sponds well to the limits we found for c 2 ia with w < — 1. 
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Figure 6. 68% and 95% confidence contours in the <^ am -w and 
c Lm"^ m pl anes using WMAP data. We see that the degeneracies 
between c? and both these parameters are very weak. The col- 
ored background indicates the mean likelihoods of the samples. 
In this model a v i s = 0. 
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Figure 5. Marginalized probability distribution for log(a„i s ) for 
a model with w = —1.2 and c? = 1. 
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4.2 Constraining dark energy sound speed 

In this section we investigate what constraints can be put 
in the dark energy sound speed, cf am , from WMAP data. 
Throughout this subsection a V i S = 0. 

From Figure [2] we see that the effect of varying cf am 
away from cf am = 1 is negligible for models with cf am > 1, 
but notable when varying cf am in the regime cf am < 1. Thus 
we would expect to find a lower limit, but not an upper limit 
on the value of cf am , and this is indeed the case. 

In Figure [5] we have plotted the 2-dimensional confi- 
dence contours in the cf am — w and cf am — f2 m planes that 
we get from the WMAP data and for a model where cf am 
is allowed to vary freely in the interval {—20, 20}. In both 
cases we see that we do get a lower limit, and also that there 
is hardly any degeneracy between cf am and these parameters 



(or any of the other parameters in our model) . Therefore, as 
in the case of c\ ia and a V i S , we would not expect to improve 
our constraints on cf am by using other types of cosmological 
data. 

Note that the intuitively plausible models with c\ am = 
reside at the lower end of a long "ridge" of c\ am values that 
fit the data almost equally well. Therefore the interpretation 
of the confidence contours should be done with some care. 
This is the reason why we do not state any lower limit on 
c fam here. Also note that the shaded backgrounds in Fig- 
ure [S] refer to the mean likelihood of the samples, whereas 
the contours show the marginalized likelihoods. For such 
highly non-gaussian distribution that we have here, there 
is a significant misma tch between these. See Appendix C in 
l|Lewis fc Bridle!l2002r i for a discussion on this. 

In Figure [7] we show the marginalized probability dis- 
tribution of cf am for the same model. In the figure we have 
also indicated the mean likelihoods of the samples by the 
background shading color. Here we clearly see the effect of 
the converging power spectrum for cf am > 1, as the likeli- 
hood for models with cf am > 1 does not decrease. We also 
see that we are allowed to have models with cf am slightly 
negative, but that the models rapidly become incompatible 
with the data for cf am < 0. 
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Figure 7. The marginalized probability distribution for c^ a 
using WMAP data. No upper limit is obtained. In this model 
ctvis = and w is allowed to vary freely. 



4.3 Effect on other parameters 

Another interesting issue is whether the extra freedom in the 
dark energy fluid will affect the constraints on the other pa- 
rameters in our cosmological model. That is, are the param- 
eter constraints in the ACDM model r obust to changes in 
(?„i„ a nd cf am . This was also studied in (|lchiki fc Takahashil 
120071 ). where they found that a V i a did not change the con- 
straints in the other cosmological parameters significantly, 
but that varying cf am shifted the other parameter ranges 
slightly. 

In Figure[8]we have plotted the marginalized likelihoods 
for different cosmological parameters in the case of a 7 pa- 
rameter model with free w but with a V i S — and cf am = 1. 
We compare this to a model where a V i 3 is allowed to vary 
freely in the interval {—20, 20} (the same model as shown in 
Figure|4]). Also shown is a model with a V i S = and cf am = 0. 
We see that the extra freedom in the a„; s parameter does 
not change the other parameter distributions significantly. 
We do however get a slight shift in the parameter distribu- 
tions by changin g cf nm from 1 to 0. This is consistent with 
the results from iflchiki fc Takahashill2007M . 

The most notable effect of changing from a model with 
= 1 to a model with cj am = 0, is that the probability 
distribution for w becomes narrower in the latter case. For 
a model with cf am = 1 we find w = {-1.47, -0.57} at 95% 
CL, while this range changes to w = { — 1.26, —0.52} for a 
model with cf om = 0. For the other parameters, the effect of 
changing cf am is not very notable. 
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5 CONSTRAINTS FROM FUTURE DATA 



As we have seen, only weak constraints can be found on 
the c^ is , a V i S and cf am parameters using present data. We 
have also argued that the effect of including other types of 
data sets, like LSS and SNIa, would not be very helpful, 
as such kinds of data sets are not affected significantly by 
these parameters. Also, unless w is significantly below -1, 
they will not serve to break parameter degeneracies for the 



Figure 8. Marginalized parameter distributions in three models 
with free w. The solid black lines correspond to a model with 



and cf 



1. 



i. The dotted red lines show the results 

lam 

from a model where a v i s is allowed to vary freely between —20 
and 20 and c? =1. The dashed blue lines show the constraints 

lam 

from a model with a v i a = and c; am = 0. 



parameters studied here. Will it then be possible to improve 
our constraints with future CMB data? 

To answer this question we have simulated a "perfect" 
CMB temperature data set, where the error bars are defined 
only from cosmic variance (CV) around a power spectrum 
generated from the best-fit ACDM model (with w = — 1, 
(? ms = and cf am = 1) from WMAP data. The likelihood 
part of CosmoMC has been modified to use this perfect data 
instead of the WMAP measurem ents. The likelihood is cal- 
culated as in l| Verde et al.ll200&t ): 



-2 log £(»»; pi) = ^(21+1) 



In 



C7I 



(7) 



For this mock data set we have used multipoles from I — 2 
to / = 2000 in our analysis, and also here we have added the 
same prior on Ho and age as earlier. A completely noise- free 
data set is of course not realistic. However, it is an inter- 
esting case to study, since effects that cannot be seen here, 
will never be possible to see using a real CMB temperature 
experiment. Note that, when using this mock data set, we 
do not include any polarization data in our analysis. 

In Figure [5] we show the constraints in the a v i s -w plane 
in a model with cf am = 1 and oi V i S < 0. This is compared 
with the results from using WMAP data alone, as also shown 
in Figure U As we can see, even in this idealized case, we do 
not see any major improvement in our constraints on a V i s . 
In this case the lower limit increases from a„,: s > —0.23 
(from WMAP data) to a vis > -0.22. 

In Figure \W\ we have used the CV limited mock data 
to redo one the most interesting case from the analysis with 
WMAP data, namely the constraints in the ci is -w plane 
with w < — 1 and <? via < 0, as shown in Figure [3] We see 
that the constraints in this area improve slightly, but not 
very significantly. Using the CV-limited mock data we find 
c vis > —17.5, compared to (? via > 19.5 using WMAP data. 

We have also redone the constraints on the dark energy 
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sound speed, cf am , using the "perfect" data set. In Figure 
llll we show the confidence contours in the cf am -w and cf am - 
O m planes, as done earlier with WMAP data. In Figure [121 
we have the corresponding 1-dimensional probability distri- 
bution. As we can see, using the perfect CMB temperature 
data, does not improve our constraints on this parameter at 
all. 

That improved CMB data does not help to constrain 
dark matter properties, is not very surprising. As seen in 
Figures Q] and effects occur mainly on large scales, where 
the error bars in the WMAP data already are dominated by 
cosmic variance. 
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Figure 11. 68% and 95% confidence contours in the cf am -w 
and c 2 am -Q m planes using only CV limited CMB temperature 
data (solid red lines). The black, dotted contours are the results 
obtained from using WMAP data, as shown in Figure \E\ The 
constraints are not improved when using "perfect" data. In this 
model a v i s = 0. 



6 DISCUSSION AND CONCLUSIONS 

We have investigated the possibility of using cosmological 
observations to constrain a possible dark energy anisotropic 
stress. We have parameterized the dark energy component 
in a very general way. Such generalized cosmological fluid is 
characterized by an equation of state (w), a speed of sound 
(cfam) an d an anisotropic stress (a V i B and cj ia ). We have 
subjected this parameterization to the most detailed and 
most extensive scrutiny this far. 

We have found that it is difficult to constrain the prop- 
erties of an imperfect dark energy fluid. If the Universe re- 
sides in the phantom regime, with w < — 1, one can put a 
weak lower limit on c^ i3 . In this case we found c^ is > —19.5 
at 95% CL. However, (? via is not bounded from above. We 
also found that in a model with negative a„„, this pa- 
rameter can be constrained to a„; s > —0.23 at 95% CL 
using WMAP data. We found no upper limit on a V i S in 
models with a V i S > 0. These results reflect the fact that a 



-1.8L 



8 Mota, Kristiansen, Koivisto, Groeneboom 




-5 5 2 10 15 20 

C lam 



Figure 12. The marginalized probability distribution for cf 
using only CV limited CMB temperature data (solid red line). 
The black, dotted line indicates the result from using WMAP 
data. The constraints are not improved by using "perfect" data. 
In this model a v i s = and w is allowed to vary freely. 



dark energy fluid with nonzero a v i s generates shear stress 
which tends to smoothen its distribution. However, the con- 
sequences to phantom dark energy are qualitatively different 
and for such a fluid, with w < —1, a shear stress drives the 
clustering. 

A similar conclusion applies to the dark energy sound 
speed, cf am , where one also can put some lower limits, but 
no upper limits on its value. 

Considering a simulated "perfect" CMB data set, only 
limited by cosmic variance, we found that improved CMB 
data sets are not likely to improve these constraints any 
further in the future. 

These results stem from the fact that there is always an 
area in the parameter space describing the possible prop- 
erties of dark energy, where the fluid is different from the 
cosmological constant or quintessence, but observationally 
indistinguishable. Optimistically though, our Universe could 
happen to reside in some other area of this parameter space. 
Indeed, in some special cases, one might be able to find a 
detect a possible imperfectness of the dark energy fluid. 

One could consider dark matter and dark energy into 
a single imperfect fluid (|Koivisto fc Motall2006l ). Since then 
the effects of shear viscosity are present already in the ear- 
lier evolution of universe, they in general have observable 
consequences. This is somewhat rel ated to models featuring 
the presence of early dark energy dDoran fc Robbers! [2009 ) 
(which is excluded by our single-parameter description of the 
evolution of the dark energy equation of state). It could be 
therefore interesting to study how the previous constraints 
derived on the presence of early dark energy will be modified 
when one takes into account possible imperfect properties of 
the fluid. 

The anisotropic stress of dark energy can have clear 
signatures in the case that it does not average exactly to zero 
at large scales. Then the nonzero a does not only appear as a 
statistically isotropic perturbation of the fluid, but will drive 
the overall expansion of the universe anisotropically. Clearly 



the resulting modificat ions of the Friedmann e quation can 
be tightly constrained l|Koivisto fc Motall2~007al '). 

However, such specific cases aside, the arguably simplest 
dark energy extension of the ACDM model introduces an 
equation of state w ^ 1. Inevitably, the more general model 
then generates dark energy perturbations, thus forcing one 
to specify also the two parameters cf am and (? vis . Again, 
arguably the simplest parameterization then considers these 
parameters constant in time. We have found then that a 
nonzero a, which is expected for any realistic fluid, typically 
escapes detection in this three-parameter model. It seems 
that it is enough to know w, and that cf am and a V i S do not 
matter as long as they are positive. 

The cosmological implication is that a canonical scalar 
field or perfect fluid representation might ultimately turn 
out to be a phenomenologically sufficient description of all 
the observational consequences of dark energy. The funda- 
mental lesson is however that even then one cannot exclude, 
by appealing to observational evidence alone, the possibility 
of imperfectness in dark energy. 
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